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Abstract 

We have investigated, using both a theoretical and an empirical approach, the frequency of low redshift 
galaxy-galaxy lensing systems in which the signature of weak lensing might be directly detectable. We find 
good agreement between these two approaches. In order to make a theoretical estimate of the weak lensing 
shear, 7 , for each galaxy in a catalogue, we have made an estimate of the asymptotic circular velocity from the 
stellar mass using three different approaches: from a simulation based relation, from an empirically-derived 
relation, and using the baryonic Tully-Fisher relation. Using data from the Galaxy and Mass Assembly 
redshift survey we estimate the frequency of detectable weak lensing at low redshift. We find that to a 
redshift of ~ 0.6, the probability of a galaxy being weakly lensed by at least 7 = 0.02 is ~ 0.01. A scatter 
in the M* — Mh relation results in a shift towards higher measured shears for a given population of galaxies. 
Given this, and the good probability of weak lensing at low redshifts, we have investigated the feasibility 
of measuring the scatter in the M* — My relation using shear statistics. This is a novel measurement, 
and is made possible because DSM is able to make individual direct shear measurements, in contrast to 
traditional weak lensing techniques which can only make statistical measurements. We estimate that for 
a shear measurement error of A 7 = 0.02 (consistent with the sensitivity of DSM), a sample of ~50,000 
spatially and spectrally resolved galaxies would allow a measurement of the scatter in the M* — Mh relation 
to be made. While there are no currently existing IFU surveys of this size, there are upcoming surveys which 
will provide this data (e.g The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), surveys with 
Hector, and the Square Kilometre Array (SKA)). 

Keywords: gravitational lensing - weak lensing 


1 Introduction 


observables. 


Weak gravitational lensing is a powerful probe of dark 
matter in the universe (eg. [Kaiser fc Squires 1993). 


Following initial investigations by Blain (2002) and 
Morales (2006), de Burgh-Day et al. (2015) have devel¬ 


oped a new method to measure the weak lensing signal 
in individual galaxies called Direct Shear Mapping 
(DSM). The primary scientific application considered 
by de Burgh-Day et al. (2015) is the measurement 
of mass and mass distribution in dark matter halos 
around individual low-redshift galaxies. In particular, 
since the dark matter halo properties can be measured 
for individual galaxies, DSM will enable the measure¬ 
ment of the dispersion in the galaxy luminous matter 
to dark matter ratio, as a function of other galaxy 
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The possibility of measuring individual galaxy dark 
matter halo masses through DSM is an exciting 
prospect, however the measurement itself is challeng¬ 
ing, and potentially observationally expensive. We have 
consequently developed the approach presented here 
for identifying the most robust candidates for such 
a measurement. We have also used this approach to 
investigate the possibility of measuring the scatter in 
the — Mh relation using shear statistics, and we 
estimate the size of the statistical sample that would 
be required to make this measurement. 

In this paper the probability of weak lensing shear has 
been estimated as a function of the redshifts of the 
source and lensing galaxies, and a catalogue of candi¬ 
date galaxy pairs is selected from the Galaxy and Mass 
Assembly Phase 1 Survey (GAMA I) Data Release 2 
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(DR2) catalogue (Driver et al. 2011 Liske et al. in 


prep). We also find that the distribution of shears in a 
galaxy sample is influenced by the relationship between 
stellar mass and halo circular velocity, and the scatter 
in this relation. With enough shear measurements it 
may be possible to constrain this relationship, and to 
measure the scatter. 


DSM uses spatially resolved velocity field information 
for an object to obtain a shear measurement from 
the velocity map. DSM assumes intrinsic rotational 
symmetry in the the velocity map, and searches for de¬ 
partures from this symmetry. This requires either radio 
data cubes or spatially resolved optical spectroscopy. 
To identify prospective targets, it is desirable to first 
obtain an estimate of the shear signal present in a 
galaxy. 


While galaxy-galaxy lensing has been used to measure 
halo masses in the past, those studies stack many 
galaxy-galaxy pairs statistically, to obtain average 
halo masses ( Brainerd et al.||1996 Hudson et al.|p!9M 


Wilson et al. 2001 Mandelbaum et al.||2006 ). In addi¬ 

tion, measurements of galaxy halo shapes have been 
made from stacked galaxy-galaxy weak lensing mea¬ 
surements ( Hoekstra et al.|2004 van Uitert et al.|2012 ). 


2 Weak lensing 


In this section we will outline the relevant weak lensing 
theory. Gravitational lensing is the deflection of light 
from some source on its path to the observer by an 
intervening mass. The deflection angle is given by 



C2 14 


( 1 ) 


where G is the gravitational constant, S is the projected 
mass distribution of the deflector, c is the speed of light 
and ^ is the distance from the deflector (i.e the impact 
parameter). Since the angle of incidence of the light has 
been altered, the source will appear to be in a different 
location to its true position. The true position of the 
source can be found by solving the lens equation 


I3 = e- ^a, (2) 

where 0 is the apparent angular separation of the de¬ 
flector and source, (3 is the true angular separation of 
the deflector and source, is the angular diameter 
distance between the deflector and source and is 
the angular diameter distance between the observer and 
source. The angular coordinates /3 and 9 can be related 
to the corresponding physical coordinates in the source 
and lens planes as 


To test our target selection algorithm, the weak 
lensing statistics of a sample of galaxies in the Galaxy 
and Mass Assembly Data Release 2 (GAMA-DR2) 
catalogue were investigated, using the stellar mass 


estimates from Taylor et al. (2011). The purpose of 


our lensing frequency algorithm is to estimate of the 
distribution of shear signals present in a dataset. 
This algorithm enables novel measurements to be 
made, and will improve the success rate of any survey 
intended to measure weak lensing via the DSM method. 

The rest of this paper is organised as follows. In Sec¬ 
tion relevant weak lensing theory is described. In 
Section a theoretical estimate of the probability of 
weak lensing at low redshift is made, following |Mort 


lock & Webster ( 2QQQ| ). In Section]^ the inputs, struc¬ 
ture and outputs of the lensing frequency algorithm are 
outlined. In Section results of the application of the 
algorithm to a dataset obtained from GAMA-DR2 are 
presented, along with an investigation of the possibility 
of using shear statistics to measure the scatter in the 
— Mh relation. Conclusions and a summary are pre¬ 
sented in Section Throughout the paper we assume a 
flat Concordance cosmology, with = 0.3, Da = 0.7, 
Hq = 70km s“^ Mpc“^, and hjo = Hq/IOO = 0.7. 


V = Dsf3 (3) 

I = D^O (4) 

(5) 

In the case of a circularly symmetric lens, and perfect 
alignment between the observer, lens and source, the 
lens equation can be solved to obtain the Einstein ra¬ 
dius, a characteristic length scale 

2 47rGM(< 6>EDd) Dds 

= -3-DA' 

where M(< ^E^d) is the mass enclosed within the 
Einstein radius, and is the observer-deflector 

angular diameter distance. 


In the weak lensing regime the light from the source 
passes well outside the Einstein radius, and the source 
is singly-imaged. In this case one can assume that so 
long as the length scale of the source is much less than 
that of the deflector, then the lensing will be linear. 
Thus it can be represented by a first-order Taylor ex¬ 
pansion, allowing Equation © to be re-expressed as a 
linear coordinate mapping between the lensed and un- 
lensed coordinate systems 


I3 = A0, 

where A is the Jacobian of transformation; 

-72 \ 

dd V -72 1 - K + 7i / ' 


(7) 

(8) 
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Here k is the convergence, and 71 and 72 are the 
two components of the shear vector 7 = ( 71 , 72 )- Equa¬ 
tion ^ can be inverted to obtain 

K = 1 — {All + ^ 22)72 ( 9 ) 

71 = -{Aii-A22)/2 ( 10 ) 

72 = —A21 = —Ai 2 . (11) 


a population of halos housing elliptical galaxies, as in 


Mortlock & Webster (2000), and a Press-Schechter 


(Press & Schechter 1974) population of halos. In both 
cases we assume an SIS halo. We will discuss the 
general steps for obtaining the expression for the weak 
lensing optical depth, and will then discuss the ellip¬ 
tical galaxy and Press-Schechter halo cases individually. 


The shear vector 7 can be rewritten in polar co¬ 
ordinates as a function of a shear magnitude, 7 , and 
an angle, 0 

7 i = 7 cos 0 ; 72 = 7 sin 0 , ( 12 ) 

where (j) is the angle of the shear vector and 
7 = \/ 7 i + 72 • work we are interested in 

the value of 7 for a given lens-source system, which is 
a function of the lens projected mass density and the 
lens-source angular separation. 


We assume a Singular Isothermal Sphere (SIS) lens pro¬ 
file throughout this paper, whose projected surface den¬ 
sity has the form 


(13) 

where is the halo velocity dispersion. Although the 
SIS profile is a primitive lens model, it is sufficient for 
this work and allows for a simple shear estimation. The 
Einstein radius for an SIS lens is given by 


oi = 




'ds 




(14) 


and from e.g. Lasky & Eluke (2009), the shear compo¬ 
nents for an SIS are 


7i 


-Dd^E 1^2 

2|^|3 


72 


-Dd^E 


^ 1 ^ 2 , 


(15) 


SO that 


D(i0e 


(16) 


Thus we have an expression for the shear magnitude as 
a function of lens-source projected separation and lens 
velocity dispersion, the latter of which can be related 
to the halo mass. 


3 An estimation of the probability of lensing 

In this section we describe the process by which a 
theoretical estimate may be made of the probability 
of a given source being weakly lensed, as a function of 
the source redshift. We begin by introducing a similar 
calculation for strong lensing in analytically solvable 
cosmologies from Mortlock & Webster (2000). We then 
adapt this work to the weak lensing case, and for a 
more realistic Concordance cosmology. In the following 
sections, we will assume two different lens populations: 


Given an estimate for the spatial distribution and mass 
of lensing galaxies in a given volume, it is possible to 
make an estimate of the distribution of shears in the 
volume. This can in turn be used to make an estimate 
of the probability distribution of weak shears across 
the sky, assuming the distribution of lensing galaxies is 
isotropic. Mortlock & Webster (2000) have used these 
arguments to make an estimate of the probability of 
lensing of quasars by elliptical galaxies for three simpli¬ 
fied cosmologies. They define the lensing optical depth, 
r, as the fraction of the source plane within which the 
lens equation has multiple solutions. It can be used as 
an estimator for strong lensing probability. The contri¬ 
bution to the total optical depth by one lensing galaxy 
is 


^/?crit 

47r 


(17) 


i.e. the fraction of the sky covered by its lensing 
cross-section. Here /3crit is the angular distance from 
the deflector where a background source transitions 
between being singly or multiply imaged. 


In the case of weak lensing, rather than being interested 
in the region where the source is multiply imaged, we 
are interested in the region where the source is singly 
imaged, yet sheared sufficiently such that it is still a 
measurable effect. This region takes the form of an an¬ 
nulus about the lens (assuming a spherically symmetric 
lens), the inner bound of which is the Einstein radius, 
6 >e, and the outer bound of which is a function of some 
limiting shear value, yum- Thus, r is re-written as a 
function of a new area, a( 7 iini), defined as the area cov¬ 
ered by this annulus: 


^g( 7 lim) — 


^(7lim) 

47r 


and 


= TT 


27ii] 


n 2 


— 1 ) 6 >e 


and substituting in equation (14) this becomes 

2 


'rg(7lim) = 


1 


27iii 


( 


f£±\ 


c) \DJ ’ 


(18) 


(19) 


( 20 ) 


where a is the velocity dispersion of the deflector, c is 
the speed of light, I^ds is the deflector-source distance, 
and Ds is the observer-source distance. 
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If it is assumed that the population of leasing objects 
are non-evolving and have uniform volume density at all 
redshifts (which is reasonable at low redshift), then the 
differential number of objects at redshift z with velocity 
dispersion a is 




dVb dUf^ 
dZ(j^ dcr 


( 21 ) 


d^^ddcr 

Here dVb /dz is the comoving volume element at redshift 


dz 


c {l + zfDA{zf 


Ho 


Eiz) 


dn, 


( 22 ) 


where c is the speed of light, Da{z) is the angular di¬ 
ameter distance at redshift z, and 


E{z) — \/(l + 
assuming a flat universe. 


(23) 


The optical depth to a redshift is then obtained by 
integrating over the optical depths of the entire popu¬ 
lation of deflectors up to Zg 

d^A^d 

T(2s,7lim)= / / -^Tg(7liin)dcrd2:d- (24) 

Jo Jo didder 

We will now discuss using this relation to estimate the 
weak leasing optical depth as a function of source red¬ 
shift for a population of elliptical galaxies, rig, and a 
population of dark matter halos, n^. 


3.1 A population of halos housing elliptical 
galaxies 


Mortlock & Webster (2000) gives the local comoving 


number density of elliptical galaxies as 
(5(l+a)-l 

exp 


drig 

* 

1 


d(T|| 

cr* ' 

Ver*/ 

where 

(T|| is 

the 


(25) 


observed line-of-sight veloc¬ 
ity dispersion, a = —1.07 ± 0.05, = (0.0019 ± 

0.003)/i7o Mpc“^, cr* = 225 ± 20 kms“^ and 
= 3.7 ± 1. The values of a and are drawn 


from Efstathiou et al. (1988), and the values of and 


S from de Vaucouleurs & Olson (1982). Substituting 


this into equation (24), the inner integral can be solved 


analytically, to give 

'^(^s:7lim) = Cg 

Jo 


(1 -h Zg)‘^ / Dd^ds 


E{zg) \ Ds 


dz. 


'g’ 


(26) 


where 

Co- = TT^ri: 


.(?) 


27iii 


-1) r(i + « + -). ( 27 ) 


For a given limiting shear value, yum, to obtain the 
probability of measuring a shear of at least the limit¬ 


ing value, one must solve equation ( [^ numerically as 
a function of the source redshift, 2 (g. 

The dashed grey line in Figure shows the lensing opti¬ 
cal depth as a function of source redshift for a limiting 
shear of yum = 0.02 based on a population of elliptical 
galaxies. The coloured lines in Figureare discussed in 
the next section. 


3.2 A Press-Schechter halo population 

The Press-Schechter mass function is given by 
dn 


__ /2 ^ dlnq-(M) ^crit 

dM “ V TTiW" dM cr(M) 


cr%t 

2cr2(M) 


(28) 


and gives the co-moving number density of dark matter 


halos as a function of halo mass (Press & Schechter 


1974). 


Substituting this into equation ( [M| ), 
^ d^rid 


T(2s,7li 


■> 77 ' 


'rd(7iim)dMd2:d, (29) 


d^ddM 

where Mmin is the mass of the smallest halo capable of 
containing a galaxy, and 

2 //7«/r 7 -1\0.316\ ^ \ 2 


7-g(7lim) = TT 7? 


1 


27iii 


- 1 


(30) 

where rj = 6.14656 x 10 ^ and we have used the relation 

(31) 


yc = 2.8xl0“2 (Mhfo^) 


from Klypin et al. ( 2011| ). Vc is the halo asymptotic 
circular velocity, and we have assumed Vc — cr. The 


lower bound of the integral over M in equation (29) has 
been truncated at M^nin because we are only interested 
in halos large enough to contain at least one galaxy. 
The value of Mmin at z 0 is usually taken to be 
Mq) 10 ( Barkana fc Loeb]|2QQl ). 


In this case, the integral over mass is not analytically 
solvable, however the two integrals are still separable, 
giving the following expression: 


(^s, 7lim) =Ch r ^ (M hfo^) dM 


f 


(1 + ( DdDds 


E{zh) V A 


(32) 


dzh, 


where 


2„2 


Ch = TT n 


27iii 


-1 - 


(33) 


and cr* ~ 321 kms“^ is a halo characteristic velocity 
dispersion, corresponding to a characteristic halo mass 
M* defined as where a{M*) = 5c{z), where 5c{z) is 
the critical density and 5c{z = 0) ~ 1.686. At z = 0, 
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Figure 1. The lensing optical depth as a function of source 
redshift and a limiting shear of 7iim = 0.02. The coloured solid 
lines show the optical depths obtained when using a population 
of lenses drawn from a Press-Schechter halo mass function. The 
three lines show the effect of different minimum halo masses. The 
minimum halo mass is usually taken to be IO^^Mq (e.g. in the 
Millennium simulation). The dashed grey line shows the optical 
depth obtained when using a population of lenses drawn from 
an elliptical galaxy population. It is reassuring to see that the 
Press-Schechter curve which best matches the dashed curve is that 
which uses the commonly used minimum halo mass of IO^^Mq. 


logio(M*/MQ) 13 ( Barkana fc LQeb|[2QQl ). 


To obtain the probability of a shear of at least yum from 
halos, as a function of source redshift, the mass and 


redshift integrals in equation (32) must be solved nu¬ 


merically. In order to illustrate the effect of varying the 


minimum halo mass, the mass integral in equation (32) 
has been solved for three values of the minimum halo 
mass: log]^o(^min/^ 0 ) = 9.5,10, and 10.5. The results 
of this are shown by the three coloured curves in 
Figure [2 As expected, the probability curve which best 
matches that of an elliptical population is for the case 
log;Lo(^min/Ar 0 ) = 10. Not Surprisingly, the higher 
the minimum halo mass the lower the probability of 
lensing, since the majority of halos are of lower mass. 


From Figure roughly 1 in 1,000 sources at z 0.2 
will be sheared by at least 7 = 0.02. In the calculation 
of halo lensing probability, four major assumptions have 
been made: 


• The lensing cross-sections of each galaxy in the pop¬ 
ulation do not overlap. This assumption breaks down 
if: 1) The redshift is high, since the number of lenses 
contributing to the optical depth increases with in¬ 
creasing source redshift. 2) Mmin is small, since this 
leads to more lenses and an increasing cross-section. 

PASA (2015) 
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3) 7 iim is small, since the cross-section of each lens is 
larger for a smaller limiting shear. 

• Both assumed halo populations do not evolve with 
redshift. This assumption is only correct over small 
redshift ranges. The elliptical and halo populations 
used in Figure are for 2 : = 0. 

• The Press-Schechter halo population is dependent on 
the minimum halo mass, i.e. the smallest halo which 
will house a galaxy. This is discussed in detail in Sec¬ 
tion |5T2 

• For this calculation we have assumed a uniform distri¬ 
bution of lens redshifts, which is not the case in real 
data. If the total lens mass distribution is specified 
however, any inhomogeneity in the galaxy distribu¬ 
tion will average out over a large sample of source 
galaxies. Thus when comparing the results of our 
theoretical calculation to lensing probabilities in real 
data, it is acceptable to use the redshifts of the source 
galaxies in the sample and treat them as homogenous. 
Such a comparison is discussed in greater detail in 
Section 15.1.11 


In spite of minor limitations, the results shown in Fig¬ 
ure are strong motivation for a more through calcula¬ 
tion and estimation of the probability of detecting weak 
lensing, since we are interested primarily in relatively 
low redshift weak lensing, where the DSM technique 
will be useful. 


4 The lensing frequency algorithm 

We now describe an algorithm for estimating the 
lensing signal per galaxy in a catalogue. We can use 
this algorithm to estimate the probability of any given 
object in the catalogue being lensed by some value, 
and to identify potentially suitable targets for lensing 
studies. 

There are three factors which will determine the 
strength of the shear signal a lens imposes on a source: 
the mass of the lens, the angular separation of the lens 
and source, and the redshifts of the lens and source. 
Different combinations of these can lead to systems in 
which the lens appears large, and almost obscuring the 
source, but has a shear of the same magnitude as well 
separated systems with a small but dense lens. The 
simplest approach to estimating the strength of the 
shear present in a potential lens-source system is to 
simply inspect their separations and redshifts. Clearly 
this is prone to misidentifications, since systems may 
be falsely rejected using this method if the lens is 
particularly dense, and systems may be falsely selected 
if the lens appears large and close to a background 
source, but is of a very low density. 
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In the lensing frequency algorithm the shear probability 
is estimated utilising information about the stellar mass 
of the lens, and the redshifts and angular separations of 
each lens-source pair. There are a number of steps in¬ 
volved in this process, beginning with the stellar masses 
of the galaxies in the catalogue, assuming a Singular 
Isothermal Sphere density distribution for each galaxy, 
estimating the halo circular velocities of each galaxy, 
and ultimately returning an estimate of the shear each 
object imposes on its nearest neighbours. These mea¬ 
surements can then be binned in shear to obtain an es¬ 
timate of the probability distribution function for shear, 
and to flag particularly promising targets for follow up 
and direct measurement of the shear. In this investi¬ 
gation we focus on galaxy-galaxy lensing only, however 
this process is equally applicable to galaxy-group lens¬ 
ing, provided good group masses are available. 


4.1 Estimation of the circular velocity, Vc 

The most uncertain step in the estimation of the esti¬ 
mated shears for each object is the calculation of the 
asymptotic circular velocity of the lens, Vc. This step 
is important as it takes us from the stellar mass to the 
total mass of each lens. Three approaches to this calcu¬ 
lation were used: 


1 . Use a power law relation between halo mass and halo 
circular velocity obtained from the Bolshoi simula¬ 
tions (Klypin et al. |MTt . In this case it is necessary 
to compute the halo masses from the stellar masses, 
which is done using the relation derived by [Moster 


et al. ( 2010 ). 


2 . Assume the circular velocity at the outermost regions 
of the disk is approximately equal to the circular ve¬ 
locity in the halo. In this case, the baryonic Tully- 
Fisher relation (McGaugh et al. 2000) can be used 
to obtain Vc from the known stellar masses of the 
objects. 

3. Use an empirically derived relation between Vc and 
(Jo from Courteau et al. (2007). This relation has a 
significant scatter. 


There are significant uncertainties in all of the meth¬ 
ods discussed here, and the discrepancies in the results 
obtained are at times significant (in particular for very 
high and very low mass objects). However, there is no 
single best method, and each method approaches the 
problem from a different starting point. These three 
methods may bracket the reality, and so by utilising all 
three methods simultaneously a good representation 
of possible values is obtained. Therefore the shear 
resulting from all three methods are computed. 

The steps taken by the algorithm are as follows: 

PASA (2015) 
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1. The objects are sorted by redshift and the angular 
size distance for each object is computed by integrat¬ 
ing over the Friedman equation: 

D = ^J[n^{l + zf + nA]-^dz. (34) 

2. Beginning with the lowest redshift galaxy, assume it 
is a lens, and perform the following steps: 


a. Compute the angular separation between the ob¬ 
ject and its neighbours out to a specified pro¬ 
jected radius ^max, excluding those at a lower 
redshift: 

^sep = yA[i?cosD]2 + A[a]2, (35) 

where A[i?cosT>] = {Ri — R 2 ) cos(T)i A D 2 ), 
AD = Di — D 2 and Ri, Di are an objects right 
ascension and declinatiorQ • 

b. Sort the neighbours by angular distance from the 
object; 

c. Compute the halo mass from the stellar mass 
( [Moster et al.||2QlQ ) 


'rn{Mh) _^( rn\ 

" vmJo 


Mh 


Ml 


+ 


mJ 


n -1 


(36) 

where is the halo mass, m is the stellar mass, 
{m/Mh)o is a normalisation. Mi is a characteris¬ 
tic mass where m{Mh)/Mh is equal to {m/Mh)o^ 
and (/) and 77 are two slopes which indicate the be¬ 
haviour of the relation at the low and high mass 
ends. The values for the free parameters used in 
this work are the best-fit values from IMoster et al.l 


( 2010 ): 


log (Ml) — 11.884_q q 23 

= 0 . 02820 t°;““i 

(/. = 1 . 0571^46 

3 + 0 . 010 . 


(37) 


V = 0 . 556 +^:^ 004 , 

d. Compute the circular velocity three different 
ways: 

i. Using a power law relation derived from 
the Bolshoi simulation ( Klypin et al.|[ 2 MT ), 
which uses a ACDM cosmology: 


Vc = 2.8 X 10“^M' 


0.316 


(38) 


where Mvir is the virial mass and it is as¬ 
sumed that Mvir = Mh. 


^The purpose of only including galaxies with ^sep < ^max in this 
step is to improve computation time, otherwise every object 
would be compared to every other higher redshift object. Back¬ 
ground objects with a large separation from the source are likely 
to have negligible shears, and so can be safely excluded. 
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ii. From the stellar mass using the baryonic 
Tully-Fisher relation ( McGaugh et al.|2QQQ ) 




0.25 


(39) 


26.25 hjQ Mq 


km ^s^. 


where A 

iii. Using an empirically derived relation be¬ 
tween Vc and (Jo (Courteau et ah ]|M7| ): 


Vc = V^O-Q. 


(40) 


These will henceforth be referred to as the 
ACDM, bTF and <Jo methods respectively, 
e. The Einstein radius of the object is computed for 
each of the three Vc values (where it is assumed 
that due to the virial theorem one can write Vc ^ 
ay), and the shear for each Vc value is computed; 

3. Move to the next highest redshift galaxy, assume it 
is a lens, and repeat the above steps. 

4. Continue iterating, increasing in redshift until the 
second-highest redshift has been reached. 


5 Application to the GAMA survey data 
release 2 catalogue 

The target selection procedure outlined in Section 
has been applied to the Galaxy and Mass Assembly 
Phase 1 Survey (GAMA I) Data Release 2 (DR 2 ) 
catalogue (Driver et al. IS Liske et al. in prep). 
The GAMA Survey is part of a larger project aiming 
to exploit the latest generation of ground-based and 
space-borne survey facilities to study cosmology and 


galaxy formation and evolution (Driver et al. 2009). 
Phase I of the GAMA Survey is a magnitude limited 
spectroscopic survey measuring galaxy spectra and 
redshifts in three equatorial regions centred at 9^, 12 ^ 
and 14.5^ (called G09, G 12 and G15 respectively), 
each with an area of 12 x 4 deg^ (Baldry et al. ]| 2010 | | . 
The fields were observed to a limiting r-band apparent 
magnitude of Tapp < 19.4, Tapp < 19.8 and Tapp < 19.4 
mag respectively. The target galaxies are distributed 
over a redshift range 0 < ^ < 0.5 with a median red¬ 
shift of 2 ) 0.17. The GAMA DR 2 catalogue contains 

all GAMA I main survey objects down to r < 19.0 
mag (for GOO and G 12 ) and r < 19.4 mag (for G15) 
including spectral redshifts (Baldry et al. 2014| ). The 
catalogue contains a total of 72,225 objects, of which 
71,599 have derived stellar masses ( Taylor et al.||20lf|^ 
GAMA survey data is available on the GAMA websit^ 

The lensing frequency algorithm has been applied to the 
GAMA I DR 2 catalogue in two ways: 


^http://www. gama-survey.org 
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1 . As the entire catalogue (with some minor cuts de¬ 
tailed below) to investigate the probability of any 
object being sheared by at least some value, and 
to identify conceivably suitable targets for potential 
follow up observation and shear measurement with 
DSM; 

2 . As a smaller sample to investigate the possibility of 

measuring the scatter in the relation, and 

to determine what sample size of galaxies is required 
to make this measurement. 

These two analyses are described in detail in the follow¬ 
ing subsections. 


5.1 Identifying weak lensing candidates in 
the GAMA survey Catalogue 

In this section we look at the entire GAMA 1 DR 2 
catalogue (with minor cuts described in the next para¬ 
graph) , to investigate the probability of detecting shear 
of at least 7 for any given object in the catalogue, and 
identify candidate targets for follow up observation and 
shear measurement with DSM. 

After removing objects with undefined or uncertain 
redshift, and selecting for stellar mass in the range 
8 < logio(M,,/M 0 ) < 12, a dataset of 69,434 objects 
was obtained. The data was passed to the lensing fre¬ 
quency algorithm in two ways; 


As the original set of galaxies (i.e. 69,434 objects). 
The shears present in the dataset were calculated via 
the three halo circular velocity estimates (AGDM, 
bTF and ao). The three — Vc relations are shown 
in Figure as dotted (ACDM), dashed (bTF) and 
solid (cTo) lines. A histogram of shears present in the 
GAMA DR 2 sample is shown in Figure as solid 
lines. The probability of a given galaxy being sheared 
by a particular value is obtained by normalising by 
the total number of galaxies in the sample. 

As a larger (synthetic) population of galaxies, to ac¬ 
count for the scatter in the — Vc relation. The 
synthetic population is produced by making 100 reali¬ 
sations of each input galaxy, with a scatter introduced 
to the realisations (giving a total dataset of 7,029,800 
objects). The scatter in the population arises from 
the intrinsic scatter in each — Vc relation and so 
is different for each of the three relations used. A scat¬ 
ter of ctacdm = 0.15 dex was introduced into Equa¬ 
tion (36) for the ACDM method ( Moster et al.| 201 Q ); 
a scatter of cTbxF = 0.14 dex was introduced into 


Equation (39) for the bTE method (McGaugh et al. 


2000 ); and a scatter of a^j^ = 0.08 dex was introduced 


into Equation (40) for the a^ method. The (Jq val¬ 
ues and associated scatter used in the ao method 
were obtained by fitting for the stellar masses us¬ 
ing velocity dispersions obtained from the Sloan Dig- 
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Figure 2. The M^c — Vc relation for three methods of calcula¬ 
tion. ‘K’ denotes the ACDM met hod ([Klypin et al. ] [2QTT1 >, ‘M’ 
the bTF metho d ([McGaugh et al.||2000| ) and ‘C’ the ao method 
( [Courteau et al.||2007| >. The dotted, dashed and solid lines show 
the relation with no scatter introduced. The blue, maroon and 
green points show a synthetic dataset obtained by making 100 
realisations of the original dataset, and introducing a scatter of 
<^ACDM = 0-15 dex, (JbTF = 0.14 dex, and a^Q = 0.08 dex for 
the ACDM ( [Moster et al.poTol , bTF ( [McGaugh et al.|2000| > and 
ctq methods respectively. The three methods agree well in the 
range 8 < log]^Q(M*/MQ) < 12 (to within ~ 15%), where most 
galaxies are situated. 



Figure 3. Histogram of estimated shears for galaxies in the 
GAMA survey. The solid, coarsely-binned lines correspond to the 
shears present in the original GAMA DR2 sample. The dotted 
lines correspond to the shears present in a synthetic dataset, ob¬ 
tained by producing 100 relisations of the original dataset, and 
introducing a scatter in the M* — Vc relation. For every galaxy in 
the GAMA survey, there are 100 galaxies in the synthetic dataset. 
Therefore for ease of comparison of the histograms of the syn¬ 
thetic and real datasets, the number of sheared objects in the 
synthetic dataset have been divided by 100. One can see that the 
three methods agree well within each of the real and synthetic 
datasets, with the synthetic datasets sitting slightly above the 
real datasets. 


ital Sky Survey Data Release 10 (SDSSDRIO; |Ahn 


|et al. [ 2 M 4 ). The — Vc relation of the synthetic 
population of galaxies is shown in Figure as blue 
(ACDM), maroon (bTF) and green (ao) points. One 
can see that in the range 8 < log]^o(^*/^o) < 
where most galaxies are situated, the three meth¬ 
ods agree well. The shears present in the synthetic 
dataset were computed for every point, and the re¬ 
sulting number of sheared objects, normalised by 100, 
are presented in Figure as dotted lines. The syn¬ 
thetic data is normalised by 100 to allow for better 
comparison to the real data, since for every galaxy 
in the real dataset there are 100 galaxies in the syn¬ 
thetic dataset. One can see that the three methods 
agree well within within each of the real and synthetic 
datasets, with the synthetic datasets sitting slightly 
above the real datasets. The implications of this up¬ 
ward shift are discussed in Section |5.2| The proba¬ 
bility of a given galaxy being sheared by a particular 
value is obtained by normalising by the total number 
of galaxies in the sample. 


5.1.1 Diseussion 

The number of galaxies with shears above the cutoff 
value of 0.02 varied between the three methods of 
calculating Vc, with the ACDM method giving the 
largest estimate of measurable shears in the data, 
and the bTF method giving the smallest estimate. 
The probability of a shear of at least 7 = 0.02 was 
P(> 7) = 0.018,0.005 and 0.007 for the ACDM, bTF 
and (Jo methods respectively. The bTF method was 
chosen for selecting candidate targets because it gives 
the smallest shear probability, and so will be the least 
likely to overestimate the shears present in a sample. 
Using the bTF method, the number of objects in 
GAMA with an estimated shear of 7 > 0.02 is 393. The 
bTF method always estimates a smaller shear than the 
ACDM and do methods for any given object, so any of 
the 393 objects selected by the bTF method would also 
have been selected by the other two methods (although 
of course the converse is not true). 

These objects were extracted from the sample, and 
matched to objects in the SDSS DRIO catalogue. The 
lens-source pairs in the sample were ranked by each 
of the authors by eye, based on their morphology, 
surface brightness, inclination angle, separation and 
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Figure 4. Two example galaxy pairs from the GAMA DR2 
Sample identified with the target selection algorithm. In both 
cases the galaxy in the crosshairs is the source galaxy (i.e. the 
galaxy being lensed). The left-hand images show thumbnails of 
the galaxy pairs from the SDSS DRIO Finding Chart Tool, while 
the middle and right-hand images show the J-band images from 
the UKIDSS survey, with the residuals from 2 D Sersic fits, taken 
from GAMA’s online Single Object Viewer tool. The top pair is at 
RA = 213.705 deg, DEC = 1.623 deg, and has an estimated shear 
of 7 = 0.023. The lens and source redshifts are 2 : = 0.128 and 
2 ; = 0.186 respectively. The bottom pair is at RA = 213.705 deg, 
DEC = 1.623 deg, and has an estimated shear of 7 = 0.053. The 
lens and source redshifts are 2 ; = 0.088 and 2 ; = 0.190 respectively. 


environment. The motivation for selecting on this 
criterion stems from the requirements of the DSM 
algorithm; that the source galaxy be undisturbed 
and stably-rotating, neither face-on nor edge-on and 
bright enough for observation with an IFU. For a more 
detailed discussion of the requirements of the DSM 
algorithm, see de Burgh-Day et ah (2015). Two highly- 
ranked example lens-source pairs are shown in Figure 


Figure shows the distribution of galaxies with 
significant estimated shear, with lens stellar mass 
plotted as a function of lens redshift. As might be 
expected, the selection favours more massive lensing 
galaxies, with significant numbers of candidates at all 
redshifts. The quality of the shear measurements how¬ 
ever depends on the characteristics of the source galaxy. 


The normalised cumulative sum was taken of the data 
shown in Figure the results of which are shown 
in Figure giving the probability of measuring a 
shear of at least 7iim- One can again see the upward 
shift in the distributions for synthetic data, relative 
to those for the real data. Also shown are two dashed 
lines corresponding to our theoretical estimate of the 
probability of lensing from Section The probability 
of at least a given shear being observed in a galaxy 
increases as the minimum shear decreases. Since the 
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logio(M.) 


Figure 5. Lens stellar mass as a function of lens redshift for 
the GAMA DR 2 sample, showing the distribution of shears with 
these parameters. Grey contours show the density of all galaxies 
in the sample. Blue squares show galaxies with a shear in the 
range 0.01 < 7 < 0.05, and red circles show those with a shear 
in the range 7 > 0.5. The shears were estimated using the bTF 
method. 


lensing cross-section of each galaxy goes as this 

is not surprising, as smaller minimum shears rapidly 
increase the fraction of the sky in which measurable 
lensing will occur. The method which gives the lowest 
overall probability is the bTF method. Since we expect 
to be able to measure shears as small as 7 = 0.02 
with DSM, we can expect to have at least a one 
in two hundred chance of measuring shear in a ran¬ 
domly chosen galaxy (from the bTF method), and up to 
a roughly one in fifty chance (from the ACDM method). 

It is obvious that the ACDM method results in much 
higher probabilities than the bTF and do methods. It is 
interesting to note that the ACDM method is the one 
which uses only simulation and theoretically-derived 
relations to obtain values of Vc from M^. In contrast, 
the bTF and do methods utilise empirically derived 
relationships between Vc and M^. One can easily 
identify the origin of the higher numbers of lensed 
galaxies from the ACDM method; the break in the 
— Vc relation for this method. While giving good 
agreement in the intermediate mass range, this results 
in a much larger corresponding halo circular velocity. 
These higher Vc galaxies will have a much larger 
lensing cross section, resulting in a larger number of 
lensed objects. Hence, while there are relatively few 
galaxies with logiQ{M^/MQ) > 12, they contribute 
strongly to the total lensing probability. This raises the 
question: how does this over-abundance of higher-He 
galaxies affect the results of ACDM-based cosmological 
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Figure 6. Estimated probability of measuring a shear of at 
least 7 for galaxies in the GAMA survey. The dotted lines are 
the probabilities with scatter introduced into the M* — Vc rela¬ 
tion (i.e. the synthetic dataset), and the solid lines are the prob¬ 
abilities with no scatter introduced. Again one can see that the 
synthetic datasets sit higher than the real datasets. The dashed 
lines are included for comparison purposes, and correspond to 
the probabilities as a function of 7 derived from the theoretical 
calculation of probability in Section using the redshifts of the 
sheared objects in the GAMA survey. 


simulations? Further consideration of this topic is 
beyond the scope of this paper, and is left to future 
work. 

As mentioned already, in Figure we have included two 
dashed lines corresponding to our theoretical estimate 
of the probability of lensing from Section The 
probability calculated in Section is the probability 
that a source at redshift 2;^ is lensed by at least some 
value 7iim- The probability obtained from the lensing 
frequency algorithm however is the probability that 
a single galaxy from the GAMA DR2 sample, chosen 
at random, is lensed by at least yum. In order to 
compare the two, we need to calculate the theoretical 
probability that one galaxy, chosen at random from a 
GAMA-like population of galaxies, is lensed by at least 
7iim- From Section the probability that some object 
in the GAMA DR2 sample is lensed by at least yum 
will be dependent on its redshift only (once a lens mass 
distribution has been assumed). Thus the probability 
that any galaxy in the sample has been lensed by at 
least 7iim will be the sum of the individual probabilities 
of each one having been lensed. Then the probability 
that one galaxy chosen at random has been lensed by 
at least yum will be the probability of any having been 
lensed, divided by the total number of galaxies in the 
sample. Figureshows the histogram of the redshifts of 
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Figure 7. Distribution of redshifts in the GAMA sample, after a 
redshift quality cut has been made. The mean and median of this 
distribution are 0.18 (solid white line) and 0.17 (dashed white 
line) respectively, and the range of the distribution is 0 < 2 : < 
0.65. 


the GAMA DR2 sample after a redshift quality cut has 
been made. While the distribution of redshifts in the 
GAMA sample is not uniform, it is sufficiently smoothly 
varying that the assumption made in Section (that 
the lenses are uniformly distributed in redshift) will 
suffice. Thus we have arrived at a method by which 
we can compare the results of Section and Section 

From Section!^ the theoretical probability of measur¬ 
ing a shear of at least 7 = 0.02 in any one galaxy in 
the GAMA survey is P(> 7) 0.002 and 0.003 for the 

halos housing elliptical galaxies and Press-Schechter 
halo populations respectively. Gomparing to the 
probabilities obtained from the GAMA DR2 sample 
(P(> 7) = 0.018,0.005 and 0.007 for the AGDM, bTF 
and (Jo methods respectively), we see the difference is 
a factor of ^ 6 for the AGDM method, and a factor of 
^ 2 for the bTF and do methods. This is an acceptable 
level of agreement, given the assumptions made in 
the calculations in Section The level of agreement 
between the two approaches can be seen in Figure 

The theoretical estimate of P(> yum) drawn from the 
Press-Schechter halo population is dependent on the 
value of the minimum halo mass in which a galaxy will 
form. As can be seen in Figure a smaller (larger) 
value of Mmin will result in a higher (lower) probability 
of measuring a shear of at least yum as a function of 
redshift. It can be seen in Figurethat adjusting Mmin 
will alter the level of agreement between the theoretical 
estimate and the probabilities derived from the GAMA 
DR2 sample. Each dashed grey line represents a mini- 


5.1 
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to place any meaningful constraints on the scatter 
in the relation. An alternative approach has been 
to utilise satellite galaxies to constrain the galaxy 
luminosity-halo mass relation. Early work involved 
stacking the central galaxies to obtain a statistical 
measure of the kinematics of the satellite galaxies 


(Erickson et al. 

1987 Zaritsky et al. 1993 

Zaritsky 

& White 

1994 Zaritsky et al. 

1997). However, recent 


work has avoided the need for stacking, largely by 
utilising the larger datasets offered by the 2dF Galaxy 


Redshift Survey (2dFGRS; Golless et al. 2001) and 
SDSS (McKay et al. 2002 Brainerd fc Specian| 


2003 


Prada et al. 2003 van den Bosch et al. |2004b|a| [Conroy 


et al. 2005[ 2007 More et al. 2009). These techniques 
to not utilise weak lensing, and do not attempt to 
measure the scatter in the — Mh relation. 


Figure 8. Comparison of the results from the lensing frequency 
algorithm with the theoretically derived Press-Schechter probabil¬ 
ities, with varying Mmin- The highest theoretically-derived proba¬ 
bility corresponds to a minimum halo mass of log]^Q(Mniin/^o) = 
8, while the lowest corresponds to log(Mmin/^o) = 12. The 
results of the lensing frequency algorithm are plotted as in Fig- 
ure[^ The bTF and cro methods best agree with a Press-Schechter 
halo population with a minimum halo mass in the range 9 < 

logio(Mmin/Mo) < 10. 


mum halo mass incremented by 0.5 logiQ{MQ). The line 
with the highest corresponding probability is that with 
log]^o(^min/Af0) = 8, aud that with the lowest corre¬ 
sponding probability log;Lo(^min/Ar0) = 12. The prob¬ 
abilities from the lensing frequency algorithm are as in 
Figure It is interesting to note that for the bTF and 
cTo methods, the best agreement with the theoretical 
approach is obtained for a minimum halo mass in the 
range 9 < log]^o(^min/Af0) < 10, whereas the currently 
accepted value is logiQ{M^i^/MQ) = 10. 


In this section we test whether DSM measurements 
of a population of galaxies could be used to measure 
the scatter in the relation. This is a novel 

measurement which is difficult with traditional weak 
lensing techniques, but is made possible with DSM be¬ 
cause it can measure individual shears with far greater 
accuracy around individual galaxies. We then briefly 
describe how to fit for the scatter in the M>^ — 
relation, and measure the scatter in a set of simu¬ 
lated shear datasets with a known scatter incorporated. 

As was noted in Section |5.1[ and can be seen in 
Figures and scatter in the relation 

results in a shift in the distribution of shears present in 
a population of galaxies towards larger shears. If the 
measurement error in the shear is sufficiently small, or 
the population of galaxies with known stellar masses is 
large enough, it is possible to measure the scatter in 
the — Mh relation by comparing the distribution of 
shears to those obtained from a relation with 

zero scatter. 


5.2 Measuring the scatter in the — Mh 
relation 

Since halos of given mass can have different halo 
concentrations, spin parameters and merger histories, 
we expect them to house galaxies with a range of 
masses. This manifests as a scatter in the stellar mass 
to halo mass relation. The value of the scatter in the 
— Mh relation is not well constrained, as directly 
measuring the masses of galaxies and their host halos 
is not trivial. Abundance matching techniques can be 
useful for describing the relationship between stellar 
mass and halo mass, however they cannot constrain the 


To perform this measurement with DSM, velocity 
maps from from an intermediate redshift survey would 
be required. The velocity maps can be obtained from 
several components of the galaxy, such as HI in radio 
wavelengths, or Integral Field Unit (IFU) observations 
of stellar velocities and gas emission in optical wave¬ 
lengths. IFU maps obtained from bright, Ha emitting 
galaxies are the most practical and easily obtainable in 
the immediate term however, and so we will focus on 
these galaxies in this analysis. 

We assume that the scatter in the M^ — Mh relation 
is lognormal, which is the standard form assumed in 


scatter in the M^ — Mh relation (Gonroy & Wechsler 

the literature (e.g. 

Behroozi et al. 

2010 Moster et al. 

2009). 

Leauthaud et al. (2012 

) used traditional weak 

2010 Behroozi et al. 2013 

). The mean of a lognormal 


lensing techniques to constrain the — Mh relation, 
however as with abundance matching they are unable 


distribution is a function of the size of the scatter 
in the distribution, and is given by exp(/i-1-cr^/2). 
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where a is the scatter in the distribution, and fi is 
the natural logarithm of the mean of the underlying 
normal distribution. It can be seen that the mean 
increases with increasing scatter. Thus if a lognormal 
scatter is introduced into the — Mh relation for a 
population of galaxies, the result is a larger number of 
galaxies with higher mass halos. It is this property of 
the lognormal scatter in the relation which 

leads to an increased probability of larger shears. 


To compare the shear distributions obtained from 
populations with and without scatter in the 
relation, the lensing frequency algorithm was applied 
to a selection of mock catalogues generated from a 
subsample of the complete GAMA DR2 sample. This 
subsample was obtained by performing the same cuts 
on the GAMA DR2 as in the previous section, along 
with the following additional cuts: O.I < 2: < 0.15, 
r < 17.5, and keeping only galaxies with H-a emission. 
The resulting dataset contained 2,861 galaxies. When 
computing shears in the mock catalogues, a maximum 
lens-source separation of T^max = 0.2 Mpc was used for 
step (2a) of the lensing frequency algorithm. 


A successful measurement of the scatter in the 
— Mh relation will be limited by survey sample size 
and the size of the scatter. To investigate the range in 
which these parameters would enable a measurement of 
the scatter, a set of simulated datasets of varying size 
was created, with a range of scatters in the 
relation. 


To produce the mock catalogues, larger populations of 
galaxies were generated from the subsample of GAMA 
DR2 galaxies by producing Monte-Carlo realisations 
of each galaxy. The resulting datasets contained 
N= 1,000, 15,000, 50,000 and 150,000 galaxies. The 
distribution of shears for a tight — Mh relation 
was computed for the simulated datasets. A scatter, 
was then included in the M* — Mh relation 
and the shears were computed again. This process was 
repeated for a selection of values of (Tm^-Mh between 
O.I and 0.5. This range of scatter is chosen to bracket 


the current literature values (e.g.. Yang et al. 


2009 


More et al.||2009] [Behroozi et al.|[2bl0[ Moster et al. 


2010). The resulting ‘true’ shears for each dataset 


were then ‘observed’ with a range of measurement 
errors, so that the ‘observed’ shear in the tight datasets 
contained a shear measurement error only, while the 
‘observed’ shear in the scattered datasets contained a 
shear measurement error, and additional scatter from 
the M^ — Mh relation. 


so? 3) What is the uncertainty on the measurement? 
To answer these questions, one first needs to establish 
that for a given number of galaxies, the distribution 
of shears arising from a tight — Mh relation can 
be distinguished from the distribution with scatter in 
the M^ — Mh relation. The distributions of shears with 
and without scatter in the M^ — Mh relation are shown 
in Figures nil Figure shows the observed shear 
as a function of the predicted shear for N = 1,000, 
15,000, 50,000 and 150,000, for a shear measurement 


of 


^ Tobs 


= 0.02 and — Mh relation scatter 


of (Jm^-Mh — 0-3 dex. The solid blue and red lines 
show the mean in the tight and scattered distributions 
respectively, and the dashed blue and red lines show 
Tier from the mean in the tight and scattered distribu¬ 
tions respectively. The shaded background shows the 
relative excess (red) or shortfall (blue) of galaxies with 
scatter relative to galaxies without scatter, in 2D bins 


of size (Ayobs x ^Jpred) = (0.02 x 0.002). Figure 10 


shows the observed shear as a function of the predicted 
shear for N = 50,000, shear measurement error of 
= 0.02, and scatters of (Jm^-Mh =0-1 and 0.5 


Tobs 


dex. Lines and shading are as in Figure Figure [TT] 
shows the observed shear for = 0.005, O.OI, 0.03, 
and O.I for N = 50,000 and o-m^-Mh — 0.3dex. Lines 
and shading are as in Figure 


The distributions of observed shears with and without 
scatter in the — Mh relation were compared using 
the Two Sample Kolmogorov-Smirnov (2SKS) Test. For 
it to be possible for the scatter in the — Mh relation 
to be measurable from comparing the distribution of 
observed to predicted shears with and without scatter, 
we require the distributions to fail the hypothesis that 
they are drawn from the same distribution, under the 
2SKS Test. That is, we require the p-value of the test 
to be small. The resulting 2SKS scores and p-values 
(denoted 5'2SKS and P2 SKs) for the combinations of n 
and (Tm^-Mh considered are presented in Table 

The next step is to fit for the scatter in — Mh rela¬ 
tion in simulated datasets. To do this, we assume that 
the observed shears are distributed according to both 
the scatter in the M^ — Mh relation, and an observa¬ 
tion error, so that the likelihood of an observed shear 
given the true underlying shear is given by 


-^(7obs|7true5 5 ^lohs) 

= J N{t - 7obs, 0, cr^„bs) (41) 

X M (r, 7t 

rue 5 


The key questions are: 1) Can we identify the scatter 
from the — Mh relation over the scatter from mea¬ 
surement error? 2) How many objects do we need to do 


where Y is a Gaussian distribution representing the 
measurement error in the DSM method, and M is a 
lognormal distribution representing the scatter in the 
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Figure 9. Distributions of observed shears as a function of predicted shears with and without scatter in the M* — relation, for 
N = 1,000, 15,000, 50,000 and 150,000 respectively (clockwise from top left), with ~ The solid blue 

and red lines show the mean in the tight and scattered distributions respectively. The dashed blue and red lines show Tier from the 
mean in the tight and scattered distributions respectively. The shaded background (and colourbar at the top of each plot) shows the 
relative excess (red) or shortfall (blue) of galaxies with scatter relative to galaxies without scatter, ri^catt ~ bins of size 

(Ayobs X Aypred) = (0.02 x 0.002). 


— Mh relation. N is given by 


by 


N{x,iJ,,a) = 


\/2p( 


TTCr 


exp 


1 f X — 11 

2 V cr 


(42) 


where /i is the mean of the distribution and a is the 
standard deviation. The lognormal, M, is a distribution 
whose logarithm is normally distributed. It is given by 


M(x, /i, a) = 


1 


\/^( 


irax 


exp 


1 (\og{x/jl) 


(43) 


where /i and a are the mean and standard deviation of 
the underlying normal distribution. The total likelihood 
of a given value of ctm^-Mh fixed then given 


k^{cFM^-Mh) = -^(7obs,i|7true,i, CF ' 

(44) 

We have fitted for the scatter in the relation 

for an assumed true scatter of (Jm^-Mh =0-3 by 
maximising log[>C(crM^-M^)] foi* a range of values 
of N and cr^^bs investigate the behaviour of the 
uncertainty in the fit with these parameters. A Fisher 
Matrix analysis was used to estimate the standard 
error in the maximum likelihood for each cr^^bs 


N. Figure 12 shows the standard error, AaM^-Mh^ 
a function of these parameters. We find that for an 
assumed scatter of (Tm^-Mu =0.3, to obtain a robust 
fit with a measurement error of cr^^bs =0.02 a dataset 
of A" ^ 50,000 DSM measurements is required, a result 
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Figure 10. Distributions of observed shears as a function of predicted shears with and without scatter in the M* — relation, for 
=0.1 dex (left plot), and 0.5 dex (right plot), with ~ ^ ~ 15,000. All lines and shaded regions are as in 

Figure 


Table 1 The Two Sample Kolmogorov-Smirnov Test score, S 2 SKS: and p-value, P 2 SKS 5 for fho scattered and tight datasets, for each 
combination of A, and cfM^—My^ investigated, n is the number of simulated datapoints for each real galaxy, N is the total number 

of simulated datapoints, and cfM^-My^ is the scatter in the M* — relation. 


N 

^Tobs 


*^28X8 

^28K8 

1,000 

0.02 

0.3 

2.5 X 10-2 

0.91 

15,000 

0.02 

0.3 

8.7 X 10-3 

0.61 

50,000 

0.02 

0.3 

1.1 X 10-2 

4.4 X 10-3 

150,000 

0.02 

0.3 

7.9 X 10-3 

1.9 X 10-3 

50,000 

0.005 

0.3 

1.6 X 10-2 

6.2 X 10-3 

50,000 

0.01 

0.3 

1.4 X 10-2 

1.4 X 10-4 

50,000 

0.03 

0.3 

1.0 X 10-2 

9.5 X 10-2 

50,000 

0.1 

0.3 

8.7 X 10-3 

4.7 X 10-2 

50,000 

0.02 

0.1 

7.6 X 10-3 

0.11 

50,000 

0.02 

0.5 

1.9 X 10-2 

1.6 X 10-3 


which is consistent with the results of the 2SKS test. 


5.2.1 Diseussion 

For fixed f^e p-value for N = 

15,000 is too large to rule out the hypothesis that the 
two samples of data come from the same distribution, 
while the p-value for N = 50,000 is sufficient to rule 


out this hypothesis, and the p-value for N = 150,000 
can easily do so. Similarly, for fixed N and (Tm^-Mu^ 
the p-value for = 0.03 is too large to confidently 
rule out the hypothesis that the two samples of data 
come from the same distribution, while the p-value for 
cr^obs ~ sufficient to rule out this hypothesis, 

and the p-value for cr^^bs =0-01 can easily do so. Not 
surprisingly, a larger scatter in the relation 
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Excess scattered galaxies, 

-150 -120 -90 -60 -30 0 30 60 90 120 150 -80 -60 -40 -20 0 20 40 60 80 



0.004 0.008 0.012 0.016 0.004 0.008 0.012 0.016 

Predicted shear, 7pped 


Figure 11. Distributions of observed shears as a function of predicted shears with and without scatter in the M* — relation, for 
cr-^obs ~ 0.005, 0.01, 0.03 and 0.1, with (Tm^-M^ — 0-3 dex and N = 15,000. All lines and shaded regions are as in Figure 


results in a smaller p-value, since a larger scatter 
directly increases the difference between the datasets 
with and without scatter in them. From the above 
considerations we conclude that in order to measure a 
scatter of (Tm^-Mh = 0-3, a sample of > 50,000 galaxies 
with a shear measurement error of ^ would 
be required. For larger values of the scatter, larger 
values of cr..y^bs smaller values of N would be 

sufficient. 


This result is confirmed in our fits for a scatter in 
the — Mh relation of 0.3 for a range of values 
of N and cr..y^bs5 presented in Figure 12 For a true 
— Mh scatter of 0.3, a shear measurement er¬ 
ror of cr^obs = 0.02, and a sample of N =50,000 shear 
measurements, we recover a fitted scatter of 0.308±0.02. 


DSM is expected to achieve measurement errors of 
cr^obs 0.02, and so shear measurement error is not 


seen to be a limiting factor in measuring scatter in the 
— Mh relation. While a sample of ^ 50, 000 galax¬ 
ies with spatial and spectral resolution does not yet 
exist, there are several surveys beginning in the near 
to intermediate future which will provide datasets of 
a sufficient size to perform this experiment, for ex¬ 
ample surveys with the Hector instrument ( Lawrenc^ 
et al.|[2M^ on the Anglo-Australian Telescope (AAT), 


and The Hobby-Eberly Telescope Dark Energy Exper¬ 
iment (HETDEX; [Hill et al.[ [2008 ), or surveys on the 
Square Kilometre Array (SKA). We conclude that while 
it would not be possible to utilise DSM to measure the 
scatter in the — Mh relation with existing lEU sur¬ 
vey data, it will be possible with data from upcoming 
surveys. 
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Figure 12. The standard error in the maximum value of 
log [>C—M^)] cLS a function of the shear measurement error, 
(j^obs values of N in the range 1,000< N <150,000. 


would be needed to measure the scatter, for a measure¬ 
ment error on the shear of 0.02 (a value consistent with 
the shear measurement accuracy achievable with DSM). 
We attempt to fit for the scatter in the rela¬ 

tion for a set of simulated datasets. The result of this 
is shown in Figure For a true scatter of 0.3, a shear 
measurement error of 0.02 and ^50,000 shear measure¬ 
ments, we recover a scatter of 0.308±0.02. It should be 
noted that the technique we have demonstrated here 
is based on a relatively untested algorithm, however it 
is adequate as an illustrative example and our results 
are promising for future measurements. While there are 
no existing IFU survey catalogues of a sufficient size 
to apply this technique, there are several surveys be¬ 
ginning in the near to intermediate future which will 
provide datasets of a sufficient size to perform this ex¬ 
periment, for example HETDEX which aims to observe 
^ 10^ galaxies, or surveys on the SKA, such as the ‘bil¬ 
lion galaxy survey’ which aims to observe ^ 10^ galax¬ 
ies. 


6 Conclusions 


We have made an analytical estimate of the frequency 
of a source being weakly lensed given a uniformly 
distributed population of lenses, following [Mortlock fc| 
Webster (2000). We have adapted their work for the 


weak leasing case, in which we consider the probability 
of the source being lensed by at least some limiting 
value 7iim. The results of this analysis suggest the 
probability of detecting weak leasing greater than a 
limiting value of yum = 0.02 in a realistically observable 
redshift range (z < 1) is non-negligible. Given this, 
we have created a leasing frequency algorithm which 
searches an input dataset for all lens-source pairs with 
an estimated shear greater than a limiting value of 
7iim- Our algorithm has been applied to a dataset 
extracted from the GAMA survey catalogue, and 
the number of objects with an estimated shear of at 
least 7 = 0.02 in the sample was found to be ^ 393. 
These targets can be matched to objects in the SDSS 
DRIO Gatalogue, and a subsample of good targets can 
be chosen from this selection for follow up observations. 


A scatter in the relation results in a shift to¬ 

wards higher measured shears for a given population of 
galaxies. Given this, we have investigated the feasibility 
of measuring the scatter in the relation using 

shear statistics. We find that for a given shear measure¬ 
ment error, our ability to differentiate between a distri¬ 
bution of shears from a ‘tight’ relation, and 

one with scatter, is dependent on the size of the scatter, 
the number of objects in the sample, and the shear mea¬ 
surement error. Eor a scatter of 0.3 dex in the 
relation, we find that a sample size of ^50,000 galaxies 
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